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Abstract 

A document has recently been posted on the arXiv [1], describing analytical formulas and results 
of particle tracking simulations, for precision tests of numerical integration algorithms for an EDM 
(electric dipole moment) storage ring. In the context of an all-electric storage ring, the authors 
cite theoretical formulas by Orlov [2], However, the reference to Orlov is to a talk at a workshop 
in 2012, and is unpublished and difficult for independent researchers to access and validate. This 
note rederives and generalizes some of Orlov’s principal results, using a Hamiltonian formalism, 
and also corrects some details in both Orlov’s note [2] and the arXiv post [1]. 

PACS numbers: 29.20.D-, 02.30.1k, 02.60.Lj 

Keywords: electric and magnetic moments, electrostatic storage rings, integrable systems, Hamiltonian 

dynamics 


* srmane001@gmail.com 


1 



A document has recently been posted on the arXiv [1], describing analytical formulas 
and results of particle tracking simulations in both magnetic and electric storage rings. The 
goal is to publish benchmarking formulas for use as precision tests of numerical integration 
algorithms for an EDM (electric dipole moment) storage ring. In the context of an all-electric 
storage ring, the authors cite theoretical formulas by Orlov [2], However, that reference is to 
an unpublished talk by Orlov at a workshop in 2012, and is difficult for readers to access and 
validate, and is moreover not a peer-reviewed reference. This note rederives and generalizes 
some of Orlov’s principal results. Admittedly this note is also not a peer-reviewed document, 
but it is publicly accessible and the contents can be independently validated. I also correct 
some details in both Orlov’s note [2] and the arXiv post [1]. 

The document [1] presents results for models of both magnetic and electric storage rings. 
Only all-electric models will be treated below. I treat a particle of mass m and charge e, 
with velocity v = (3c and Lorentz factor 7 = l/^/l — /3 2 , moving in a prescribed external 
electrostatic field E. The speed of light is set to unity below c = 1. The specific model of 
an all-electric ring treated in [1] is a homogenous weak focusing ring. (See Section 5 of [1].) 
I employ cylindrical coordinates (r, 9, z ) and the design radius is denoted by r 0 . (Similarly 
the reference values of other quantities are also denoted by a subscript “0” e.g. 70 , etc.) 
The field index n is defined so that, in the median plane, the radial electric field component 
is E r oc 1 /r 1+n . The two cases studied in [1] are that of no vertical focusing (cylindrical 
capacitor, purely radial electric field, field index n — 0) and very weak vertical focusing 
(field index 0 < n 1). The authors in [1] denote the vertical direction by y and employ 
the notation 9 y for what I call d(z/ro)/d9 , the slope of the vertical motion. 

For the case of no vertical focusing, the orbit is a vertical spiral with a constant pitch 
angle (i.e. dz/d.9 or 9 y is constant). It is stated in Section 5.1 in [1] that “Y. Orlov[15] solved 
the orbital motion for an electrostatic field with no focusing. In this case, the estimates for 
the average values of A 7/70 and A r/r 0 take the following form:” [i.e. eqs. (23) and (24) in 


[ 1 ]] 



( 1 ) 


( 2 ) 


Note that actually no average (• • •) is required on the values of 7 and Ar in this model. 
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“Ref. 15” in the above statement is Orlov’s note [2], In fact Orlov did not derive the 
above expressions. (Curiously, Orlov is a coauthor of the arXiv post [1].) Equations (1) and 
(2) were derived by Ivan Koop, whose analysis was reproduced in an Appendix in a paper I 
published [4], with Koop’s kind permission [5]. See eqs. (A5) and (A7) in [4], respectively, for 
the above expressions. I have personally verified the correctness of Koop’s elegant solution. 

For the case of nonzero vertical focusing, the authors state (Section 5.2 in [1], the authors 
denote the field index by m) “With weak focusing such that 0 < m 1, the parameters 
analytically estimated by Y. Orlov[15] are given by Equations 25 and 26 below: [i.e. eqs. (25) 
and (26) in [1]] 

(t1> = °’ (3) 

= (4) 


which hold for times much larger than the period of vertical oscillations.” Recall Ref. 15 
in [1] is Orlov’s note [2], In fact, only eq. (3) appears in [2], The second expression eq. (4) 
was derived in a later note by Orlov [3] (also unpublished). Hence of the four formulas 
eqs. (l)-(4), only eq. (3) appears in Orlov’s note [2], 

Nevertheless, the expressions in eqs. (3) and (4) are correct and I derive (and generalize) 
them below. For contact with the analyses in [1] and [2], the independent variable is the 
time t. The coordinates are (r,9,z), the conjugate momenta are ( p r ,Po,Pz ) and I define 
x — r — r 0 and p x — p r . The model is an all-electric homogenous weak focusing ring, with 
an electrostatic potential <h(r, z). The Hamiltonian is 


H = 


2 I 2 I 2 I Pe 
m +p r +p z + — 


1/2 


+ e$(r, z). 


( 5 ) 


Then H does not depend on 6 hence pg is conserved. There is no rf cavity, so H is also an 
integral of the motion. The equations of motion are 

(6a) 
(6b) 
(6c) 
(6d) 


dx 
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For the case n — 0 (no vertical focusing), the potential is logarithmic e$ = myo/do ln(r/r 0 ). 
However, the analysis below treats bounded vertical oscillations, where the field index is 
n > 0. The potential is given by a hypergeometric function, and to the required order 


e$ = 


wyoA 


n 

m 7o/3 ( 


Tq „ /n n 1 2 " 




n 


- "i7o/?o 


1 - 




2 2 2 


n 2 z 2 

~ ^T 2+ "' 


n + lf n z 2 


( 7 ) 




2 r 2 
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72+7; j 1_ ( n + 2 )~ 


r 0 


Define cuo = Pq/tq as the angular revolution frequency. For the vertical motion, to linear 
order 


dp z 

dt 




■ nz = — Po^o n ■ 


r 0 


Also dz/dt ~ p z /Ho so 


C^2 

dt 2 


~ —cn 0 n 2 . 


( 8 ) 


( 9 ) 


It is well known that the vertical betatron tune is given by u z = n. We parameterize the 
vertical betatron oscillations using an amplitude parameter z' 0 and an initial phase 0 2O 


— = — sin(u z u 0 t + 0 2O ) , 

r 0 v z 

z 1 = — = z' 0 cos(u z u> 0 t + (j ) z0 ). 
Po 


(10a) 

(10b) 


Next we treat the horizontal motion. The model treated in [2] is that there are no free radial 
or longitudinal oscillations. The radial and longitudinal motions are driven by the coupling 
to the vertical oscillations. Normally, to linear order, we say that the horizontal and vertical 
motions are uncoupled, but in this analysis we include coupling terms of 0(z 2 ) and 0(p 2 ), 
e.g. see the expression for the potential $ in eq. (7) above. This means x = 0 (zq ) and 
Px = O(zq) are of the second order in small quantities. I also set H = H 0 (l + AH/Ho), 
where AH / Ho is also of the second order in small quantities. Then 
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Then dx/dt ~ p x /H 0 and 

dp x Pe/r 2 1 r2 
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This yields 


d 2 x 2 
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It is well known that the horizontal betatron tune is given by u 2 = 2 — /3q — n. Since the 
radial motion consists of bounded oscillations, one must have ( dp x /dt) = 0. We use the 
result (( Pz/Po ) 2 ) = u l({ z / r o) 2 ) = n{{z/r 0 ) 2 ) below. Hence 
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Equation (4) and Orlov’s result in [3] are special cases of the above for AH/Ho = 0. Note 
also that 
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It follows that 
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Equation (3) and Orlov’s result in [2] are special cases of the above for AH/Ho = 0. 
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• Orlov derived eqs. (3) and (4) in [2] and [3] respectively, but under the approximation 
of very weak vertical focusing 0 < n 1. (This is also stated in [1].) 1 found that such 
a restriction is unnecessary: the above derivation did not require any conditions on the 
held index other than n > 0 (so as to have bounded vertical oscillations). Tracking 
simulations confirm that eqs. (3) and (4) are valid for arbitrary values 0 < n < 1. 
There may be a caveat that the value of v z should not be rational, to avoid orbital 
resonances, to justify the statistical averages. 

• As stated by Orlov [2], the above results are also valid in the presence of rf and 
synchrotron oscillations. We can see this as follows: in the presence of rf, the only 
change is that the value of H is not constant but the average is (H) = H$. We just 
substitute (A H/H 0 ) = 0 in the above derivation, and the above expressions for the 
averages will follow. Tracking simulations confirm eqs. (3) and (4) are valid in the 
presence of rf. 

• As also stated by Orlov [2], the above derivations assume the motion is driven entirely 
by the vertical betatron oscillations (p. 14 in [2]): “Assume the presence of only vertical 
oscillations and the fields as in eq. (5). The energy in (2) depends only on y and, in 
accordance with this, ...” Also p. 2 in [3]: “It is taken into account that the shifts 
5x (of the radius) and S'f are quadratic or of a higher-order effect (in the presence of 
RF).” As I explained above, the coupling terms in the equations for the radial motion 
are of 0(z 2 ) and 0{jp 2 ). The same remark applies to the synchrotron oscillations. 

• It is stated after eqs. (25) and (26) in [1] (i.e. eqs. (3) and (4) above) “Note that these 
values depend only on the pitch angle, not on the ring geometry, ...” It is not clear 
that eqs. (3) and (4) do not depend on the ring geometry. Orlov stated (p. 15 in 
[2]) “Taking into account that (dy / dr) 2 /(/Ty)g = $1 = m(y/R) 2 , on the average, ...” 
Here r is the proper time and R is the design ring radius. In terms of my notation, 
this states that (z' 2 ) = n((^/r 0 ) 2 ) = $q. Such a result is valid only in a homogenous 
weak focusing ring. In general, in a model with bends and straight sections (e.g. drift 
spaces), the vertical betatron oscillations must be parameterized using Twiss functions 
and the Courant-Snyder invariant must appear in the formulas, and also v 2 ^ n. Hence 
it is not proved that eqs. (3) and (4) do not depend on the ring geometry. 
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It should be possible to employ the above formalism also for the case of no vertical 
focusing, i.e. n — 0. The answer is instructive. The solution for the vertical motion in this 
case is 

— = Z 0 uj 0 t , — = z 0 . (18) 

r 0 Po 

The potential is independent of z so we set n[z/r Q ) 2 = 0. We set AH/ Hq = 0 for now (this is 
important). Then setting dp x /dt = 0 in eq. (12) (and recall that no average (• • •) is required 
on the values of x, etc.) yields 

x 1 

r 0 


pI = j2 7o 

32 ^2 _ .2 


(19) 


2-ftjPo 7o + 1 

This is not the same as Koop’s result in eq. (2). The reason is that H ^ Ho in Koop’s 
derivation. Using eq. (1) and (2) and a logarithmic potential and ln(r/r 0 ) — A r/r 0 yields 
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It is well known that for n — 0, the centripetal condition for a circular (spiral) orbit depends 
only on the kinetic energy and not on the orbit radius. The orbit radius is determined by 
the potential energy. Koop treats an orbit which is synchronous with the reference particle, 
and the consequence is that H ^ H 0 . 

However, the value of A 7/70 for a spiral orbit does not depend on the potential energy, 
hence the above formalism should yield eq. (1). Let us verify this. Using eq. (19), 
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( 21 ) 


Ho Ho 
This agrees with eq. (1), as required. 

In Fig. 9 in [1], the authors claim close agreement of their numerical tracking simulation 
results with eq. (1) (top graph) and eq. ( 2 ) (bottom graph). Also in Section 5.2 in [ 1 ], the 
authors state “There is an apparent gap between Equations 23 and 24 and Equations 25 and 
26 in the limit as m —)■ 0. The transition between focusing and no focusing can exist since 
the latter formulas hold only for averages over times much larger than the period of vertical 
oscillations.” Here the field index is denoted by m and “Equations 23 and 24” (resp. 25 and 
26) are the formulas without (resp. with) vertical focusing. However, the situation is more 
subtle. The formulas with vertical focusing (eqs. (3) and (4)) are derived with (H) = H 0 . 
There is no cognizance in [1] that for a model with no vertical focusing, only the value of 
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A 7/70 is uniquely determined (and the expression in eq. ( 1 ) is independent of the value of 
AH/Hq), while the value of At/tq is arbitrary and is determined by the potential energy 
(and H ^ H 0 in eq. (2)). In fact for a model without vertical focusing, the exact value of 
the radius of a spiral orbit is 


L = e (H-mr,)/(H 0 ft) _ g(A/r,/ H 0 -A r ho),lft 

r 0 


( 22 ) 


Note that AH/H 0 and A 7/70 are not required to be small. 

Also in Orlov’s derivations in [2] and [3], the orbit path length L is not changed from the 
reference value (i.e. L = L$ after one turn), if the motion is driven by vertical oscillations, 
i.e. there is vertical focusing. This is not the case when there is no vertical focusing and the 
orbit is synchronous with the reference particle, which is the case in Koop’s derivation of 
eq. (2). The path length after one turn is 
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This fact does not seem to be noted in [1], nor in Orlov’s analyses [2, 3]. 


3~7o 2 
7o 2 + 1 


(23) 
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